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We consider an isotropic Fermi liquid in two dimensions near the n = 2 Pomeranchuk instability in 
the charge channel. The order parameter is a quadrupolar stress tensor with two bosonic shear modes 
with polarizations longitudinal and transverse to the quadrupolar momentum tensor. Longitudinal 
and transverse bosonic modes are characterized by dynamical exponents zy — 3 and z± = 2, 
respectively. Previous studies have found that such a system exhibits multiscale quantum criticality 
with two different energy scales to ~ ^"^11.^ , where ^ is the correlation length. We study the impact 
of the multiple energy scales on the electron Green function. The interaction with the critical z« — 3 



mode is known to give rise to a local self-energy that develops a non-Fermi liquid form, T,{uj) 



j2/3 

for frequencies larger than the energy scale lo ~ ^~'^. We find that the exchange of transverse z± — 2 
fluctuations leads to a logarithmically singular renormalizations of the quasiparticle residue Z and 
the vertex F. We derive and solve renormalization group equations for the flow of Z and F and show 
that the system develops an anomalous dimension at the nematic quantum-critical point (QCP). 
As a result, the spectral function at a fixed to and varying k has a non-Lorentzian form. Away from 
the QCP, we find that the flow of Z is cut at the energy scale ufl oc ^^^, associated with the z — 1 
dynamics of electrons. The z± = 2 energy scale, uj ~ £,~^, affects the flow of Z only if one includes 
into the theory self-interaction of transverse fluctuations. 

PACS numbers: 



I. INTRODUCTION 

The behavior of Fermi liquids (FL) near Pomeranchuk 
instabihties attracted high interest in the last few years 
because of a generic theoretical interest and potential ap- 
plications for the cupratesi and ruthenatesi^ In the FL 
notations, a Pomeranchuk instability occurs when one 
of the harmonics g^,™ of the Landau quasiparticle in- 
teraction function approaches — 1 (a = c, s stands for 
charge or spin, and n is the value of the angular momen- 
tum). Examples of Pomeranchuk instabilities include 
phase separation (gc.o = ^l)i ^ ferromagnetic transi- 
tion {gs,o = —1), at which the Fermi surfaces (FS) of 
spin- up and spin-down fermions split apart, and nematic- 
type transitions in the charge^"— and spin channels j^^""— 
which lower the rotational symmetry of the FS. 

The subject of this paper is the d— wave (n — 2) charge 
nematic instability in an isotropic Fermi liquid at T = 0, 
i.e., a charge nematic quantum-critical point (QCP). A 
charge Pomeranchuk instability in the d— wave charge 
channel was first introduced by Halboth and Metzner— 
and Yamase and Kohno^ in the context of the RG analy- 
sis of potential instabilities of a 2D Hubbard model. The 
d— wave instability in an isotropic system was analyzed 
by Oganesyan et al^ and in a number of later papersi^rJ^ 
The order parameter for a d— wave charge nematic tran- 
sition is the expectation value of the quadrupolar elec- 



momentum tensor (cos 2^ and sin2(/) terms). These two 
modes arc characterized by different dynamicsj^ The lon- 
gitudinal mode is Landau-ovcrdampcd by particle-hole 
pairs and has a dynamical exponent 2|| ~ 3. The trans- 
verse mode, on the other hand, remains undamped with 
Z[ =2. Each mode has an associated "bosonic mass- 



shell" energy scale. 
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^-"=11 and 



-'FL 
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related 



to the correlation length, ^. The nematic instability is 
thus a quantum phase transition with multiple energy 
scales. 

Multiscale criticality of the nematic transition was re- 
cently analyzed by Zacharias et al. (Ref. Il2[ ) within a 
bosonic Ginzburg-Landau </)^-theory. These authors con- 
sidered the interplay of the two modes and the resulting 
manifestations of the multiple energy scales in thermody- 
namics. Longitudinal fluctuations have larger dynamical 
exponent, zh = 3, i.e., larger phase space, and domi- 
nate the critical specific heat C{T), which undergoes a 
crossover from a FL behavior at T < lo'L to C(T) ^ T'^/^ 



for T > 



■'FL- 



Transverse fluctuations with zj^ = 2 have 



tron density c\Q 



^(^jjOq, 
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2didj, and 



smaller phase space and account only for subleading cor- 
rections to C{T). At the same time, the effective dimen- 
sion of the transverse mode d + z± = 4: is upper critical, 
and its self-interaction gives rise to singular logarithmic 
corrections to thermodynamics. As a result, critical ther- 
modynamics becomes sensitive to the second energy scale 
as well. In particular, Zacharias et al. argued that 



U! 



FL 



i,j = X, y. It can be interpreted as a traceless quadrupo- 
lar stress tensor representing elastic shear modes of the 
Fermi surface (FS) (Refs. 1 511 121). In spatial dimensions 
d = 2, there exist two bosonic shear modes with polar- 
izations longitudinal and transverse to the quadrupolar 



the temperature dependence of the correlation length ^ 
changes qualitatively at a temperature T ■ 



^FL- 



In the present work we evaluate the electron self- 
energy, I](k, Wm), close to the nematic transition focusing 
on the interaction with transverse fluctuations. The self- 



energy correction, Ey , due to the exchange of longitudinal 
bosons has been calculated earlier <^"—ii£iii^ Its most sin- 
gular part only depends on frequency and is thus purely 
local. For frequencies \uJm\ < ^fl' ^11 i*^ linear in ujm as 
for a FL liquid but obeys E||(a;m) 



2 /3 

u!m for higher fre- 



2/3 



quencies. The Wm behavior of the self-energy is directly 
related to T'^'^ behavior of the specific heat. Recent ad- 
vanced studies of longitudinal fluctuations at a nematic 
transition found that (i) there is infinite number of uim 
terms with relative coefficients 0(1), even if the theory 
is extended to large N (the planar diagrams, Ref. [20) 

2/3 

such that there is no guarantee that cum behavior sur- 
vives, and (ii) there arc singular logarithmical corrections 
at third and higher-orders of the loop expansiou i^^'^^ In 
this paper we show that another set of logarithmically 
singular corrections arises already in lowest loop order 
due to interaction with transverse fluctuations. We show 
that the latter renormalizes the residue of the fermionic 
propagator contributing to the anomalous dimension for 
the fermions. 

The issue of rcnormalizations by transverse fluctua- 
tions is tricky. As the contribution of the transverse fluc- 
tuations to C(T) is negligible except for the renormal- 
ization of the correlation length, one might expect that 
the electron self-energy remains unaffected by the inter- 
action with the transverse mode. Indeed, we find that the 
self-energy is smooth at the transverse "mass-shell" en- 
ergy scale WpL. We argue, however, that the interaction 
of electrons with transverse fluctuations brings about a 
third energy scale wfl ~ vp^~^, where vp is the Fermi 
velocity. This scale corresponds to z = 1 fermionic exci- 
tations near a fermionic mass-shell. We find that above 
this scale the exchange of transverse bosons gives rise to 
a singular logarithmic correction to the residue Z of the 
electron Green's function. Summing up the leading loga- 
rithmic contributions, we obtain that they exponentiate 
and contribute to an anomalous dimension if the distance 
to the fermionic mass-shell exceeds this third energy scale 
wfl- At criticality, wpL = and at any distance to the 
fermionic mass-shell, the electron Green function then 
acquires the form 



G(k,w,„) 



1 



isignuj„i\uj,n\'^/^ujy^ 



£k 



l-rr 



(1.1) 



where wq ~ Ep /{kpci)'^ , Ep and kp is the Fermi en- 
ergy and momentum, respectively, and a is a quadrupolar 
scattering length. The anomalous dimension, 77, at one- 
loop order is attributed to the transversal fluctuations, 
771L = J]j_, and we find rij_ = l/{2kpa). 

In order to keep our calculations under control we as- 
sumed that the interaction is sufficiently long-ranged in 
real space such that a is much larger than the inverse 
average distance between particles, i.e. kpa » 1, so that 
ri± is small. It was shown by Metlitski and Sachdev^l 
and Mross et al^ that the longitudinal fluctuations also 
contribute to the anomalous dimension, ry, but this con- 



tribution appears at third loop-order and is beyond the 
accuracy of the present study. 

The singular behavior of Z affects the spectral function 
A(k, w), particularly the momentum distribution curve 
(MDC) measured in ARPES experiments at a fixed, 
small uj and varying momentum k. In the presence of 
a singular fermionic residue Z, the momentum tails of 
the MDC arc no longer of Lorentzian, 1/e^ form, but 
rather behave as 



AiKuj) 
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at Efc > lui 
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(1.2) 



This is the experimentally detectable prediction of the 
theory. 

The input for our calculations is the assumption that 
a nematic critical point does exist, i.e., that there is no 
pre-emptive instability at some finite correlation length 
^. A pre-emptive pairing instability is always a possi- 
bility, but the corresponding Tc is generally quite low.— 
Metlitski and Sachdev pointed out^i another potential 
pre-emptive instability: a non-singular correction to the 
q^ momentum dependence of the static bosonic propa- 
gator (which is the same in longitudinal and transverse 
channels) is large in a large A'^ limit and may lead to a 
spiral- type instability already at a finite ^3^ This issue 
was further discussed by Mross et al^ who demonstrated 
that such a correction is small and under control if the 
static bosonic propagator has the form 1/q^^'^ instead of 
l/q^, and e is small (the non-trivial limit considered in 
Ref. I2I is e < 1, Af > 1, eAf = 0(1)). Wc assume in this 
paper that the regime of small e extends to e = 1, A^ = 1, 
i.e, that there is no pre-emptive spiral-type instability. 

The paper is organized as follows. In the next Section 
we introduce the effective fermion-boson model for the 
charge nematic quantum-critical point. In Section [llll we 
analyze the renormalizations arising from the exchange 
of longitudinal fluctuations and particularly focus on how 
the longitudinal mode influences the dynamics of the 
transverse fluctuations. In Section IIVI wc consider the 
renormalizations due to transverse fluctuations and de- 
rive the result for the transverse self-energy that leads 
to Eq. p.ip . The paper concludes with a summary and 
discussion of the results. 



II. FERMION-BOSON MODEL FOR THE 

CHARGE NEMATIC QUANTUM-CRITICAL 

POINT 

Effective fermion-boson models near quantum-critical 
points have been discussed in the literature for ne- 
matio^'i'^'i^ and other—— cases, and wc simply state 
the result: the proper low-energy theory near a charge 
QGP is obtained by integrating out high-energy degrees 
of freedom and is described by an effective Hamilto- 
nian with a four-fermion interaction mediated by a static 
propagator of soft bosons. This reflects the fact that the 
transition itself and the propagator of soft static bosons 



arc produced by fcrmions with high energies, of order Ep. 
The fermionic self-energy E(k, Wm) is sensitive to the dy- 
namics of bosons, which by virtue of an energy conserva- 
tion law comes from low-energy fermions and has to be 
calculated within the low-energy model, self-consistently 
with S](k, Wm)- 

Because there are two soft boson modes near a nematic 
transition in d = 2, the low-energy theory is a trace of 
a 2 X 2 matrixi^ii^ Non-diagonal terms of this matrix 
do not play a role near a QCP and we neglect them. 
The two diagonal terms describe interactions mediated 
by longitudinal and transverse bosons and are given by 

i?ll = (2.1) 

Yl Xst(q)4,q4,q4-hq/2,aCp-q/2,/3Cp+q/2,/jCk-q/2,a, 
k,p,q 

H^ = (2.2) 

/ . Xst(q)rfk,qC'p.qCk+q/2,«Cp_q/2,/3'^P+q/2,/3Ck-q/2.Q, 
k,p,q 



where 



^k,q 



%/2cos(20k,q), 4.q - %/2sin(2,^k,q) (2.3) 



are "d— wave"formfactors, and the angles, (/)k,q = Z(k, q) 
and 0p,q = Z(p,q), are between the bosonic momen- 
tum q and the two fermionic momenta, k and p, respec- 
tively. We assume, as in earlier works^iii^ that the static 
bosonic propagator Xsti^) is analytic at small q and is 
given by 



Xst(q) = 



Xo 



1 + ffc.2 + (aq)^ + 



(2.4) 



where xo is the product of the square of electron-boson 
coupling and the 2D electron density of states, v = m/7r, 
(within RPA, xo = 1/i^j Refll3l and for simplicity we 
assume below that this holds, i.e., xo^ — !)• As we 
already mentioned above, we assume that the effective 
quadrupole interaction is sufficiently long-ranged in real 
space, so that we can treat l/{kpa) as a small parameter. 
The nematic instability occurs when the Landau pa- 
rameter gc,2 in Eq. p.4|) reaches —1. Using p.4|) . we can 
introduce a correlation length 



C 



(2.5) 



v/l+5c,2 
It diverges as the transition is approached, gc.2 — > — 1. 

A. Dynamics of critical nematic fluctuations 

The polarization of electrons gives rise to a retardation 
of the effective interaction, and, as a result, the suscep- 
tibilities in Eqs. (|2.ip and (|2.2p acquire dynamic parts: 



Xa(q,^r: 



1 



1 



I' 1 + 5c,2 + (ag)2 + 6Ila (q, ^m)/!^ 



(2.6) 




ib) 




FIG. 1: Lowest order (a) transverse, H^, and (b) longitudinal, 
II|i, polarizations. 



where a =|| , _L, and the dynamic part of the polarization 
bubble is defined as SUa{q,ilm) = ^a{<i,^7n)^^a{(l,0)- 
The polarizations for free electrons are given by 



ni°nq,a: 



(2.7) 



-E Z_^ ('^k,q) Go,k-|-q/2,w„+0„/2Go,k-q/2,( 



-0„/2, 



Eb is the bare electron Green 



where Gq (k, utm) = i^ 
function, see Fig. [T] Interestingly, the dynamics dis- 
tinguishes between fluctuations that are longitudinal or 
transverse to the quadrupolar momentum tensor of the 
electrons)^ To see this, it is convenient to integrate first 
over £fe. The requirement that the poles of the two Green 
functions have to be in different half-planes for the Sk 
integration automatically restricts the integration over 
fermionic frequency, Wm , to a range of the size of the ex- 
ternal bosonic frequency flm- The dynamical SHq, then 
takes the form 



mW(q,f}™) = ^^J-„f^ 
where the functions J-a are defined as 



(2., 



-^a(s) = / 



6k,q V^^-1 



2tt s ~ cos ( 



'k,q 



Evaluating the angular integrals one obtains 



J'll(s) = -2(2s2-i) 2s 



1 - 2s2 ls + 1 



1 + ,s V s - 1 



jr^(s) = -4s 1 - 2s^ + 2s(s - 1) 



At small s, 

J^ll (s) = — 2i sign(/77is) + 4s + 
J^±{s) ~ —4s — 8i s sign(/?7is) 



s + 1 
s- 1 



(2.9) 

(2.10) 
(2.11) 



(2.12) 
(2.13) 



Both J-i\{s) and J'±{s) contain branch cuts originating 
from the pole in (|2.8p . In the longitudinal case, the 
branch-cut non-analyticity determines the leading behav- 
ior in the limit \flm\ <^ vpq: J^|| (s) ~ —2i sign(/TOs), giv- 
ing rise to a dynamical term \nm\/(vFq) in SIl.. (q, flm), 
characteristic for Landau damping. For the transversal 
polarization, on the other hand, the small frequency limit 



is analytic: ^±(s) ~ —4s, because the transverse form 
factor d^ = •\/2sin(20k.q) vanishes when both electrons 
are on mass-shell, i.e., 4>k.q ~ tt/2. This limiting behav- 
ior can be directly obtained from Eq. (|2.9p by expand- 
ing the integrand to first order in s, corresponding to 
the approximation of a quasi-static virtual particle-hole 
pair. The angular integral in Eq. (|2.9p then averages the 
direction of the ccntcr-of-mass momentum of the pair, 
k, over the FS. As a result, the J7,„ dependence of the 
transverse polarization JII^ (q, flm) starts quadratically 
in frequency, {iljn/{vF(l))'^- The branch cut of Eq. (j2.1ip 
yields a damping term in SII^^'{q, fl^) only at the third 
order in rim/ivFl)- This damping term then can be 
safely neglected in the critical scaling limit \ilm\ ^ {vrq)- 
The lowest order polarizations thus yield different dy- 
namics for the transverse and longitudinal susceptibilities 
in the limit \flm\ ^ vpq 



XO||(q,^m) 



^l+5c,2 + M2 + 2^ 



\n„ 



pq 



Xoj_(q, f^m) = - 



^t 



^i + .gc,2 + M2 + 4(^ 



(2.14) 

(2.15) 



As a consequence, the theory close to quantum critical- 
ity is characterized by multiple dynamical exponents z^ 
Whereas the non-analytic ilm dependence of the longi- 
tudinal mode corresponds to Landau damping charac- 
terized by an exponent Z|| ~ 3, the transversal mode is 
instead propagating with z± ^ 2. 



III. LONGITUDINAL FLUCTUATIONS 

When electrons scatter off critical nematic fluctua- 
tions, singular FL corrections arise. In the present sec- 
tion, we concentrate on the exchange of longitudinal fluc- 
tuations. We flrst briefly review the results for the "lon- 
gitudinal" self-energy and then discuss how singular self- 
energy and vertex corrections due to longitudinal fluc- 
tuations affect the polarizations Slia- Of our particu- 
lar interest here is the renormalization of the transverse 
(5n_L which we then use to calculate the "transverse" self- 
energy, which is the main subject of our work. 



A. Longitudinal self-energy 

The one-loop electronic self-energy arising from the ex- 
change of a longitudinal boson, see Fig. ^a.), has been 
calculated beforei^'^'i^ For completeness, we repeat the 
result here: 



E||(k,a: 



E 

q,f2™ 



(d, 



I 
k+q/2,q^ 



'Gck+q.w,, 



hn„Xo|l,q,o„ 
(3.1) 




FIG. 2: Lowest order (a) self-energy and (b) vertex correction 
due to the exchange of a longitudinal boson (dashed line). The 
curly line at the vertex represents either a longitudinal or a 
transverse boson. 



Neglecting in the form factor d" the bosonic momentum 
q compared to the fermionic k and expanding the de- 
nominator of Go in Eq. (j3.ip up to first order in q, the 
self-energy takes the form 



E||(k,w„) 



dqq 

27r 



rff^mXO||(q,^m) 



27r 



vpq 



X Th 



i{LOra +^m) - £fc 



vpq 



(3.2) 



where J-\\{s) is given by (|2.10p . The behavior of 
E||(k,aj,„) at low energies is determined by the limit- 
ing behavior of the T\\ for small argument: J^ii (s) « 
— 2i sign(/ms). Substituting this form, we obtain the 
local self-energy E||(w,„), independent of momentum k. 



S||(tJm) 



< LO- 



FL 



isigncj,„|w„p/3cJo'^ if|i:^m| > i^fl' 



(3.3) 



where ojo ~ Ep/{kpa)'^. At small frequencies, E|| has a 
Fermi-liquid form with the prefactor 



A = 



1 



e 



2{kpa)y'l 4- 5c, 2 2kpa'^ 



(3.4) 



that diverges at criticality. This FL frequency range 
shrinks, upon approaching criticality, as 



Ep 



CO- 



FL 



{kpa)^X^ 



(xC 



-3 



(3.5) 



The dependence Wp^ (x ^ "^ is characteristic for zy = 3 

the exchange of a 



> U!- 



FL' 



dynamics. For frequencies |a;„ 

2/3 

longitudinal boson leads to a non- Fermi liquid, cum form 
of the self-energy. 

The momentum-dependent part of Sii comes from 
higher terms in the expansion of J^ii (s) in s in (|3.2[) 
and is regular, T,\\(k,uim = 0) oc Cfc, and the prefactor 
is small in l/(fci^a)J^ Three-loop diagrams give rise to 
Ck log Ck terms in the self-energy, and eventually have to 
be taken into consideration . ^^'^^ For simplicity, we re- 
strict ourselves here to the one-loop longitudinal self- 
energy Eq. (j3.3|) . Within this approximation the elec- 
tron Green's function dressed by longitudinal fluctua- 



tions reads 

G||(k,w™) 



lUJr, 



S||(Wm) 



£k 



— ki?(k — lip)/'!' 



(3.6) 



where by virtue of purely local I]|[(a;,„), Zn = m/m* = 
[1 + 9E||(a;m)/9(iWm)]~"'". In the FL regime, Zy = 

m/m* = 1/(1 + A), in the non-FL regime cu > Wp^, both 
Zu and m,* /m become functions of frequency. 





(a) 



FIG. 3: Lowest order (a) self-energy and (b) vertex correction 
to the polarizations from the exchange of a longitudinal bo- 
son (dashed hne). The outgoing curly lines represent either a 
longitudinal or transverse boson for the respective Ha- Dia- 
gram (a) is redundant if the one-loop polarizations in Fig. [T] 
are evaluated with dressed Green functions (13.61). 



B. Renormalization of the polarizations 



The lowest-order longitudinal self-energy p.3p implies 
the breakdown of the FL close to quantum criticality. 
The question arises as to how this singular FL correc- 
tion feeds back into the polarizations and thus modifies 
the dynamics of the boson propagators, Eqs. (|2.14p and 

Previous studies have established that the Landau 
damping term of the longitudinal propagator is robust 
against dressing by longitudinal bosons by two reasons J^ 
First, dressing fermions in the polarization bubble by lon- 
gitudinal self-energy adds to the \Qm\/{vFq) term the 
overall factor (Znm* /m.)''^, which remains equal to one 
even when Z|| vanishes and m,* /m diverges. Second, ver- 
tex corrections need to be evaluated at k = kp, uj„i = 
^m = 0, and 9 — >■ 0, because the Landau damping term 
is the leading term in an expansion in fln/ivpq)- In this 
limit, vertex corrections are small in l/(kpa) and are thus 
irrelevant. As a result, the Landau danrping remains uir- 
modified. 

The situation is more tricky for the polarization of 
the transverse bosons. The dressing of fernrions in 
the transversal polarization IIj^ yields for the dynam- 
ical il'^/{vpq)'^ term an extra factor Z'n{m.* /m)^ = 

1 + A which diverges upon approaching the QCPJ^ At 
the QCP, the transverse bubble evaluated with dressed 
fernrions has a non-analytic \^m\^ /{vpq)'^ form. Si- 
multaneously, however, vertex corrections become rel- 
evant because in, e.g., the FL regime the leading 
vertex correction at small but finite Xflm is 6T ^ 
^m / {vpq cos (j)). The integrand of 11^ in Eq. (|2.9p also 
contains flm sin 2cf>/ (vpq cos (j>) ternr, and the product of 
the two yields il.f^/{vFq)'^ term in SII±, which has the 
same functional form as the term coming from the renor- 
malizations of the electron propagators. 

This interplay between vertex renormalizations and 
dressing up of fermions (i.e., self-energy corrections) has 
been considered perturbatively by Zacharias et al. (Ref. 
Il2f ). They explicitly computed one- loop self-energy and 
vertex correction terms with the assumption that the cur- 
vature of the FS can be neglected. They found that these 
two corrections exactly cancel each other, i.e., 5Il± pre- 
serves the same form as for free fermions. This cancel- 
lation is a somewhat unexpected result. Vertex and self- 



energy corrections do cancel each other exactly in the op- 
posite limit q = and flm — )- 0, as required by the Ward 
identity associated with the conservation of the particle 
number. However, there is no generic requirement that 
there must be a cancellation at arbitrary Q.m/(vpq), and, 
in particular, at \U,m\/vFq <C 1. Yet, a very similar can- 
cellation between self-energy and vertex corrections to 
the polarization bubble at small ^m/ivpq) has been re- 
ported by Kim et al. for the problem of fermions coupled 
to gauge fluctuations^^ 

Below we reconsider this problenr with particular em- 
phasis on the role of the curvature of the FS and higher- 
order terms. Although our primary interest is the trans- 
verse polarization, for conrpletcness we present calcula- 
tions for both transverse and longitudinal bubbles. For 
the transverse bubble, we show that singular self-energy 
and vertex corrections cancel, in agreement with Ref. Il2l . 
and the fully renormalized dynamical SIl± preserves the 
same n'^/{vpq)'^ form as the bare SIl±. Moreover, the 
the renormalization of the prefactor is small in l/{kpa). 
For longitudinal polarization, we indeed flnd that the 
leading Landau-damping term is not renormalized. The 
subleading, ^"^/{vpqY term is renormalized, but the 
renormalization is again small in l/{kpa). These results 
mean that both longitudinal and transverse fluctuations 
maintain their FL form and are robust against the dress- 
ing of polarizations by longitudinal bosons. In the follow- 
ing, we first consider the dressing of 5IVa by longitudinal 
bosons in the small Vim/ {vpq) limit. It turns out that 
in this limit, 5T\a is obtained by a systematic perturba- 
tive loop expansion. In a second step, we evaluate 5\la 
in the FL regime for an arbitrary value of fl„i/i'vpq), by 
summing up ladder series of vertex correction diagrams. 



1. Polarization dressed with longitudinal fluctuations, small 

^m/ivFq)- 

We first evaluate the polarizations in Fig. [T] with the 
renornralized electron Green's function, Eq. p.6p . Sub- 
stituting (|3.6p into (|2.7p and integrating over nronrenta 



transverse to the FS we obtain, instead of Eq. (|2.8I 
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duJrn^a 

vpq J V vpq 
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En„ 



n„/2 



duJn 



rf0k,c 



'^l^ 



-f2„,/2 



^ 2tt En^.cu^-VFqcos(j>i^^q 



where we introduced the abbreviation E< 



n„ 



, (3.7) 






^||w„+nm/2~S||LJm-n„/2j and the index R indicates that 
dressed electron Green functions, G||, are used. 

It is convenient to substract from this expression the 
bare polarization SHa evaluated with the undressed 
Green's functions. The difference is 



(3.8) 
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At small fl„i/{vFq) (or, more accurately, at small 
En,n.ij^/{vFq)), the frequency-dependent terms in the 
denominator can be neglected, and evaluating the inte- 
gral over (^k,q we obtain 



5n«it)(q,a„)-<5nW(q,an) 



(3.9) 



' {vpqf 



dlOri 



i1^l2 



^||cj„+n™/2 



-am/2j 1 



where 6|| = 1 and 6^ = —1. 

This correction can in fact be attributed to the per- 
turbative diagram presented in Fig. [31[a). Evaluating 
the frequency integral, we find that in the FL regime 
|r2m| "C WpL, the integral in (|3.9p is of order \^\^l{vFqf 
and at larger |ri„j| '^ uj\^ it is of order \^^^l'^ jivFq)^ ■ 
Whereas such a non-analytic correction is only sublead- 
ing for the longitudinal polarization, it would dominate 
over the Vt^^jiypqf' term of the transversal one in the 
limit \Vlra\ <C Vpq- However, this singular renormaliza- 
tion coming from the self-energy insertion is canceled out 
with the renormalization coming from the vertex correc- 
tion, as we now demonstrate. The polarization bubble 
with a vertex correction is presented in Fig. Eljb). Eval- 
uating this diagram, we obtain 

XVl^\qMrn)^\ J2 {dl^T ST^,^,M,^„,) (3.10) 



/? 



k,LjjTi ,a 



G'||,k-Hq/2,w„+n„/2G'||^k-q/2,Ld„-n„/2- 

where Jrk,(^„(q, r^™) is shown in Fig. [2][b). We assume 
and then verify that dV only depends on the orientation 
but not on the magnitude of fermionic momentum k, so 



that we can explicitly integrate over e^ in Eq. p.lO[) : this 
automatically restricts the oj^ integral to a range of size 
rim. We obtain 



du^^\q,n^) 
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ii' I dujj-i 

'-n™/2 



(3.11) 
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.(q,^rn) 



-. 27r En„ 
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Wc now compute (5rk,c^„(q, 5^m)- Neglecting small 
bosonic compared to large fermionic momenta in form 
factors and using Green functions dressed with the lon- 
gitudinal self-energy, we obtain 

<S(q-^"0 = ^ E (4.q')%o||(q',^:„) (3.12) 

q',0;„ 
'^ G'||^k+q'+q/2,a;„-f-n;„+0,„/2G||.k+q'-q/2,LJ,„+0;„-n„/2- 

The integration over the orientation of the bosonic mo- 
mentum q' is dominated by angles for which one of the 
two Green functions is on mass-shell. Keeping only the 
contributions attributed to these poles Eq. p.l2|) be- 
comes 
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^>F J 27r Jo 
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En^,uj„, -Vpq COS (pu^q 



27r ./o 27r 



dq' / / o' ^ 



(£^o,„,c^,„ - Vpq COS (j)k^qf - (^ sin(^k,qj 



(3.13) 



The term in the denominator involving sin 0k, q is due 
to the curvature of the FS, i.e., due to the fact that the 
dispersion ek^+q = WF9C0S(/)k,q -I- q^ /{2m). 
Expression (j3.13p can be further simplified to 

<S(q'^™) = ^'1;:"'^"'"^' : ^J':::?"'^' (3-14) 



1^!. 



En,r,,uj„, -Vpq cos (j)i^^q 
VFq\sin(f)i^^q\ 



2{kpaY Eq^^^^ - WFg cos </)k,q 
1 



En„,,u,„ - Vpq cos 0k,q + I 2(fcfa)2A ^^g'^"'^ 



with A defined in Eq. p.4p . The first term on the right 
hand side of Eq. p.l4p is the result for the vertex cor- 
rection if we neglect the curvature of the FS. It is pro- 
portional to the difference of longitudinal self-energies 
at frequencies of the intermediate fermions. Subsituting 
this term into (|3.1ip . neglecting En,„,ujm compared to 
f i?? cos ^k.q and integrating over 0k, q, we find that it ex- 
actly cancels out the renormalization of the polarization 
bubble coming from self-energy insertions, Eq. p.9p . 

The remaining term in Eq. p.l4p is attributed to the 
FS curvature. It is determined by the pole arising from 
the presence of a finite curvature term in Eq. (|3.13p . 
i.e., by putting the intermediate particle- hole pair on- 
shell. For its evaluation, one could actually neglect 






FIG. 4: Ladder approximation for the vertex summing up 
the repeated exchange of longitudinal bosons (dashed line). 
The curly line represents either a transverse or a longitudinal 
boson. 



FIG. 5: Polarization containing a vertex in the ladder approx- 
imation, see Fig. ID 



2. Polarizations at arbitrary 0,,n/{vFq): ladder 
approximation 



the frequency dependence of the longitudinal susceptibil- 
ity, xo||(q:^m)- The overall Qm dependence is just the 
phase-space to create a particle-hole pair. Substituting 
this term into p. Ill) , again negecting Eq,^,^^ compared 
to VFqcos(j)\^^f^, and integrating over (?!)k,q, we obtain the 
final result for the leading vertex correction in the limit 



<5nfl3t^^)(q,a^) = -6„ 



Av 



nl 



[kpay {vpqY' 



(3.15) 



where, we remind, 6q = 1 for the longitudinal polariza- 
tion, and 5q = — 1 for the transverse polarization. 

Note that to obtain this result one could just ne- 
glect the small imaginary part in the denominator in 
the last line of Eq. p.l4|) . The integral over (/)k,q is 
then determined by the third-order pole l/(iO signf2„j — 
WF? cos 0k,q)^ m Eq. ([3ll|) once we substitute ST^^"! 
into this formula. This physically corresponds to putting 
all four intermediate fermions on mass-shell, the asso- 
ciated phase-space constraint explains the (0,^/ (ypq))"^ 
dependence in Eq. (|3.15p . 

One can straightforwardly check that higher-order 
vertex correction terms only contribute higher-powers 
of rim/wpg. i.e., Eq. (|3.15p is the full result for the 
interaction-induced renormalization of the polarization 
bubble to order {i^m/vrq)^- For the longitudinal polar- 
ization, this renormalization is subleading compared to 
the \Qrn\/(vFq) Landau-damping term. For the trans- 
verse polarization, Eq. p.l5|) has the same fl'^^/{vpq)'^ 
dependence as the frec-fernmion bubble, Eq. (|2.15p . i.e., 
the exchange by longitudinal nematic fluctuations does 
affect the prefactor for the transverse polarization bub- 
ble. Still, the correction to the prefactor is small in 
l/(fcFa)^, and can therefore be safely neglected. 



Note that Eq. (|3.15l) is independent of the actual fre- 
quency dependence of the dressed fermionic Green's func- 
tion, the only thing that matters is that the sign of the 
dynamic part of the dressed G^^(k, Wm) is the same as 
the sign of uJm- In fact, Eq. (|3.15p could be also derived 
using a quasi-static approximation for the fermion prop- 
agator, G'~^(k, Wm) Ri iO signaj„j —Sk or using bare Green 
functions. Gq, instead of G. 



We now extend the analysis to arbitrary flm/iypq)- 
To avoid cumbersome expressions, we restrict the con- 
sideration to the FL regime, \fl,rn\ ^ '^fl- "^^^ extention 
to non-FL regime at larger fi^ is straightforward^^ but 
requires more efforts. 

In the FL regime, we can use the FL form of the lon- 
gitudinal self-energy, E||(a;m) = AiWm, and the one-loop 
vertex correction becomes 



<S'(q."™) 



(3.16) 



iSim.A 



(1 + \)i^m - w_Fgcos(/)k,q + i^^fj^^^j^j^signil. 



fori 



series 



Higer-order vertex corrections 

((5r^_^^(q,an))" and remain 0(1) when A|a«| > vpq. 
We verified that at kpa ^ 1, which we assume to hold, 
leading vertex corrections form a ladder series, while 
non-ladder terms are small in l/lkpa) (RefflSl). The 
ladder diagrams form a geometrical series, see Fig. 21 
and can be easily summed up. The resulting vertex in 
the ladder approximation reads 



^laddc 
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(3.17) 
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The polarizations are now obtained from the diagram 
with dressed fermion lines and the full vertex, see Fig. [SJ 

n„(q,a„) = i J2 Kq)'r'^'''"(k,q,an) (3.18) 
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For its dynamical part we obtain 

dlla{ci,iljn) ^ in^v (3.19) 
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^'^ iVlm - Vpq (cos (/)k,q - i 2[kfayi\ I sin (/ik,q 
iQ.m{y + A) - Vpq fcOS0k,q -«2i^^§^|si^'^k,q| 



irim(l + A) - u_Fgcos(/)k,q 



For q — and finite J7m, this reduces to the density of 
states, 6Ua = i', as it sfiould be by the Ward identity. 
At arbitrary ilm/ivpq), the full SHa can be expressed as 
the sum of two terms, 



SIla{q,^r, 



(3.20) 



6U(^\ci,n^) + -^(^) V^ 



irir, 



{kpaY \VFqJ \VFq 
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Close to the QCP A is large, and the functions V are 
given by 





(b) 



FIG. 6: Lowest order (a) self-energy and (b) vertex correction 
due to the exchange of a transverse boson (wiggly line). The 
curly line at the vertex represents either a longitudinal or a 
transverse boson. 



Vo.{s,X) 



i signlms 



dlW I sin 01 



(s — cos(/))2(s(l + A) — cos( 



(3.21) 



where, we remind, s = iflm/vpci, '^||(</') = 
d±{4') = V^sincf). At small s, i.e., at 



: v^cos 



< 



lims^o^||(s,A) 



-1 and lims_j.o ^±(sj A) = 1, 



and 

vpq, 
and 



Eq. p.20p reproduces the perturbative result p.l5|) . At 
finite s, there is an intermediate regime at large A when 
\s\ <^ 1, |s|(l + A) ~ 1, where Va{s, A) « 'Pq(sA) becomes 
a function of sA, e.g., V±{sX) oc l/^/l + (isA)^. At \s\ ^ 
1 ^ |s|A, both IIj^ and H n beh ave as Pa{s, A) ex l/(sA), 
e.g., the second term in (J3.21I) has the same \n,n\/vFq, 
form as Landau damping, bnt is suppressed by 1/A. Fi- 
nally, for large s we get Va{s, A) ex l/(s'^A). 

In the non-FL quantum-critical regime the computa- 
tion of SIla{ci,^m) becomes more complex, but a com- 
parison with the FL result shows that at S||(nm) ^ vpq 
and flm ^ vpq, 



<5nW(q,r!„o 



(3.22) 
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The overall conclusion of the analysis in this section 
is that, as long as we approximate the self-energy by 
Ti\\{uJm), the original forms of both polarizations, lon- 
gitudinal and transverse, are preserved. In particular, 
neglecting quantitative corrections that are small in the 
parameter l/(fci?a), the full Xaiq^^m) can be safely ap- 
proximated by Eqs. (|2.14p and (|2.15l) . 



IV. TRANSVERSE FLUCTUATIONS 

We now turn to the renormalizations arising from the 
exchange of critical transverse fluctuations, i.e., bosonic 
nematic modes with momentum q and a polarization 
transverse to the corresponding quadrupolar momentum 
tensor, 5ijq^ — 2qiqj. We first analyze perturbative one- 
loop corrections to the electron self-energy and the ver- 
tex and demonstrate that they are logarithmically large. 
Afterwards, we sum up the leading logarithmic singular- 
ities and show that the full result is the appearance of an 
anomalous dimension of the electron Green function. 



A. One-loop self-energy 

The lowest order electron self-energy arising from the 
exchange of a transverse boson is shown in Fig. EJa). 
Following the same steps leading to Eq. p. 21) . we obtain 
for the transverse self-energy 
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X Ti 
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dqq 
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d^7n XO±(q, ^m) 



^" J-oc- 27r 
Vpq 
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(4.1) 



where J^± is given by Eq. (j2.1ip . and the index R again 
indicates that electron Green functions dressed with Sy 
were used for its evaluation. Consider first the contribu- 
tion obtained by setting external energy and momentum 
exactly on mass-shell, iuJm + ^\\(j-^m) ~ £k = 0. In the 
low-energy limit, we approximate J\l by its value at a 
small argument and obtain 



^fliEli) 
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rf^m Xo^ (q, ^m) 
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iypq? 



!l"m-Hn„ 



The integration over q is straighforward and perform- 
ing it and then integrating over frequency, we obtain 
that one-loop transverse self-energy on the mass shell be- 
haves as WmA^/'^/(fci?a)^/'^ at the lowest energies and as 

Wm Wq /{kpa) at the QCP. (In both limits the Q,m inte- 
gral samples the frequency regime where the longitudinal 
self-energy has the non-FL form, E||(wm) ^ \^m\'^^^-) 
This additional self-energy correction is smaller than 
Sj[(a;m) both in the FL regime and at the QCP, although 
at the QCP the relative smallness is only in l/{kpa). We 
see therefore that mass-shell transverse self-energy is es- 
sentially irrelevant. In fact, we will show in Sec. IIVC 11 
below that the full transverse self-energy on the mass 



shell is even smaller as E^*^^ (k, Wm) 



mass shell 



ally canceled by other terms. 

The important transversal self-energy (|4.ip is its off- 
shell part. Expanding Eq. (|4.1[) in the distance to the 



mass-shell, (5S = E — E| 



shell 



wc obtain 



5I]^'-''^(k,w,„) = - (icjm + S||(w,„) - Efe) / — — 
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27r {vpqY 



vpq 



(4.3) 



where J^'j_ is a derivative of J^± with respect to its argu- 
ment. The frequency integral is determined by the pole of 
Xo± and the remaining momentum integral is logarithmi- 
cally large. To logarithmic accuracy, the lower boundary 
of the momentum integral is determined by the distance 
to the fermionic mass-shell Wi^gmin — |*'^m + S||('^rn)~£fc| 
such that an expansion of Eq. (|4.ip remains justified. At 
large momenta, the integral is cut by the T^ function. 
Effectively, the upper boundary is given by Qmax — 1/a 
where Qm/ivFq), taken at the bosonic pole, becomes 
of order one. In this momentum and frequency range, 
^m/{vFq) dominates over the difference of longitudinal 
self-energies divided by vpq in the argument of the J-"^ 
function, the latter being at most of order l/{kpa) <C 1. 
So we can use a small-argument expansion for J^^ with 
-^j_(0) = ~4. We then obtain, with logarithmic accuracy 
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X log 



-'""' + "^11 ^""^"^VKO) (4.4) 
8kpa 
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FIG. 7: Two loop fermionic self-energy diagrams with the ex- 
change of a transverse (wiggly line) and longitudinal boson 
(dashed line). Diagram (a) is redundant if the one- loop dia- 
gram in Fig. [S{a) is evaluated with a dressed electron Green 
function (13. 6t. 



and the transverse susceptibility instead of the longitu- 
dinal ones. The leading contribution to the transvese 
vertex correction again comes from the regime where 
rim + 5]||(r2m) ^ vpq^ and the two fermionic Green's 
functions can be approximated by their static forms 
Gk+q.Lj„ ~ l/(~WF'7COS(?!)k,q). Substituting these forms 
into the vertex correction diagram, we obtain 
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(q, ^^ 



dKn 
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dq'q' 

~2^ 



xo_L(g',f^™) 



2 sin^ 2(t> 



^ 2tt {vpq' cos (j^y 



(4.6) 



This momentum regime is bounded from be- 
low by the distance to the two fermionic 
mass shells of the internal fermions, WF9min — 
max{|iWm±ir2„,/2-|-S]||j^^^±n^/2 -£k±q/2|}- The 

upper bound gmax — l/a again ensures that the 
internal bosonic energy taken at its pole remains suf- 
ficiently small. The remaining integral in Eq. (|4.6p is 
logarithmically large. 



The transversal self-energy thus yields a logarithmically 
singular correction to the fermionic residue Z. The log- 
arithm is cut-off at the scale. 



Ep 



WFL 



{kpayX 



vpr'- 



(4.5) 



The relation wfl oc ^~^ is characteristic for the z = 1 
dynamics of electrons. The corresponding momentum 
scale (7fl = topi^/vp ex ^~^. The z = \ scaling actually 
comes from the upper limit of the logarithmical integral 
for which, as we said, V,„i/{vpq) ^ [aqY = C'(l). 



B. One-loop vertex correction 

We next show that not only the self-energy but also the 
vertex correction due to the exchange of a transverse bo- 
son is logarithmically enhanced. The vertex correction 
is presented in Fig. EKb). It is given by the expression 
similar to Eq. (|3.12p . but with transverse form factors 



<S(q,f^™) = -^x 



2kpa 



(4.7) 



Ep/{kpa) 



ma,x{\iuJm ± «^m/2 -I- ^\\uj^±n^/2 - ek±q/2|,WFL} 



log 



where cjfl is defined in Eq. (|4.5 



C. Two- loop corrections 

In the next loop-order, there are mixed diagrams that 
involve transverse and longitudinal bosons, and diagrams 
with two transverse bosons. For the self-energy, we 
demonstrate that mixed diagrams just cancel the domi- 
nant term in the on-shell one- loop self-energy, Eq. (|4.2p . 
We show that the self-energy and the vertex corrections 
containing two transverse boson exchanges contribute 
log^ terms. We also discuss the cancellation of singular 
logarithms in the polarizations. 
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1. Two-loop self-energy with transverse and longitudinal 
boson exchange 

Consider the two-loop diagrams to the self-energy in 
Fig. [7] obtained by the exchange of a single transverse 
and a single longitudinal boson. Diagram (a) can be dis- 
regarded; it is redundant as we already evaluated the 
one- loop diagram in Fig. ^a) with a dressed Green func- 
tion G||. So we need to evaluate the remaining three: 
diagram (b) is a longitudinal self-energy diagram with 
a transversal self-energy insertion, and (c) and (d) are 
self-energy diagrams with vertex corrections. 

Consider first the diagram in Fig.[7Jb). We have 

S«0')(k,f7„) = i ^ {dl^r (G||,k+q..„H-n„)' (4.8) 



/S 



iStn 



^ ^_Lk+q.w,„+s:2„Xo||,q,n„- 

Instead of using the result Eq. (|4.4|) for T.j_ here, it 
is actually more convenient to perform first the angular 
integration over the orientation of bosonic momentum q. 
The resulting expression can then be approximated as 
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The two diagrams in Fig.[71[c) and (d) are equivalent and 
can be expressed as 






flSt„ 
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where 6T 



R^) 



k+q/2,^„+n„/2(q'^™) is given by (EH. To 
leading order in l/(fci?a), we can neglect the curvature 
and approximate this vertex by the first term on the right 
hand side of (|3.14p . To the same accuracy, we can neglect 
^r2m,wm in the denominator of p.l4[) (i.e., approximate 
the propagator of an intermediate particlc-holc pair in 
the vertex correction by a static limit), and also approx- 
imate the fermion Green's function by its static limit. 
Applying these approximations, wc obtain 



dqq 1"°° dn„i XO-L(q:^ri 
27r y_o, 2tt {vFqf 



(4.11) 
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As announced, the sum of the two-loop contributions, 
Eqs. (|4.9p and (|4.1ip . cancels the leading term in the 
dressed one-loop transverse self-energy on the mass-shell, 
Eq. (|4.2p . The sublcading terms are not canceled, but 
they are small and of no relevance. 




ib) 




FIG. 8: Two loop fermionic self-energy diagrams with the 
exchange of two transverse bosons (wiggly line) . 



2. Two-loop self-energy with exchange of two transverse 
bosons 



We now turn to the two-loop diagrams for the self- 
energy containing two transverse boson exchanges, see 
Fig. [51 We verified that, to logarithmic accuracy, the 
self-energy diagram (a) is dominated by the regime where 
the momentum of the inner transverse boson is larger 
than the momentum of the outer one. Using the one- 
loop result (|4.4p for the internal self-energy, we obtain 



1 1 
2k^'p 



X;-Ri-P)(k, ri„j) - ^p;^^ 2_/ ('^k,q)^'^l|k+q,u;„+O„X0±q,O„ 

Ef /{kpa) 



X log 



max{|iw„ -Hif}„ + S||^^+n,„ -ek+ql^^^FL} 

(4.12) 



Diagram (b) has two dominant contributions: from the 
regime where the first boson is larger than the second one 
and vice versa. Both contributions can be re-expressed 
via the one-loop vertex correction due to transverse boson 
exchange. Using Eq. (14. 7p for this vertex correction we 
obtain 

I]^lSl>)(k,a„) = --^i ^ (4.q)'G||k+q..„.+O„X0±q.O„ 



X log 



qJ2„ 

Ep/{kpa) 



max {|iw„, + i^m + S||„^+n,„ - Ek+ql, ^fl} 

(4.13) 



Observe that the "vertex" diagram (b) thus gives twice 
the contribution of the "self-energy" diagram (a), see 
Eq. (|4.12p . and is of opposite sign. The two diagrams 
then partially cancel each other, and the full result is 

onehalf of S^lS')(k,a„)- 

For the remaining integration, we focus on the off- 
shell part and expand the Green's function to linear 
order in the distance to the mass-shell. The integra- 
tion over frequency Q,m is dominated by the bosonic 
pole of the transverse susceptibility. The leading con- 
tribution to the momentum integral comes from the re- 
gion where the momentum component, q_\_^ perpendicu- 
lar to the external fermionic momentum is much larger 
than £,~^ ■ In this regime, the integrand depends on the 
longitudinal component, gii , only via the form factor, 
(dj^ )^ = '&q\q\/{q\ + 9m)^j and the integration over q\\^ 
is easily performed. The final integration over q± then 
reduces to j dq±lq±\ogq^^. CoUccing the prefactors, we 
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FIG. 10: Correction to the polarizations due to the exchange 
of a transverse boson (wiggly hne). The external curly line 
represents either a longitudinal or transverse boson. 



FIG. 9: Two-loop vertex corrections with exchange of two 
transverse bosons (wiggly line). 



then obtain a log correction to the residue Z, 



5E«lSt^'')(k,an) 



->m + S|tw„ 
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2(2fcj.a)2 



(4.14) 



ma.x{\iuJrn + ^\\{uJm) - ek\,UJpi^} 



3. Two-loop vertex correction 

We next consider logarithmic singularities in the vertex 
correction at two loop order, arising from the exchange 
of transverse bosons. The corresponding diagrams arc 
presented in Fig. [9l We find partial cancellations be- 
tween the contributions, similar to the cancellations be- 
tween the two-loop diagrams for the self-energy. Using 
one-loop results for the diagrams in Fig. [6l we obtained 
that the diagrams (a) and (c) with self-energy insertions 
cancel the diagrams (b) and (d) with vertex correction in- 
sertions. Out of two remaining diagrams, the important 
one is the diagram (e) with the two ladder-type vertex 
corrections. For the inner vertex correction we can use 
the one- loop result (|4.7p , and the remaining integrals are 
evaluated in a manner similar to how Eq. (j4.14p was 
obtained. Performing the calculation, we obtain 



<.S(q,^™) 



log' 



1 



2{2kFaf 

Ep/{kpa) 



(4.15) 



max {\iuj,n ± ^ + S||a;,„±0„/2 



^k±q/2 



I.'^fl} 



The remaining diagram (f) does not have a log term 
- it vanishes after the integration over the direction of 
bosonic momenta. As a result, Eq. (|4.15p is the full result 
for the vertex correction at two-loop order. 



4-. Polarization corrections 

Because the exchange of transverse bosons leads to 
logarithmic singularities in the electron self-energy, the 
question arises whether it affects similarly strongly elec- 
tron polarizations. If it does, logarithmical renormaliza- 
tions of the self-energies and the polarizations have to be 



analyzed self-consistently and in parallel. We show, how- 
ever, that this is not the case. In the evaluation of the 
self-energy we have already witnessed a partial cancella- 
tion of logarithmic singularities arising from self-energy 
insertions and vertex corrections. We demonstrate below 
that in the polarization bubble, self-energy and vertex 
corrections due to a transverse boson exchange conspire 
so that singular logarithmic terms are exactly canceled 
out. 

The polarization bubble with transverse self-energy 
and vertex corrections is shown in Fig. 1101 Express- 
ing each of the two diagrams in terms of fcrmionic and 
bosonic propagators and performing calculations in the 
same way as before we find that each diagram is log- 
arithmically divergent. However, comparing them, we 
find after simple manipulations that the integrands are 
identical and only differ in the overall sign. The issue is 
then the interplay between the combinatorial factors for 
the two diagrams. The self-energy can be inserted into 
the upper or the lower fermion propagator, hence the di- 
agram (a) has a combinatorial factor of two. The vertex 
correction diagram (b) does not contain this factor, how- 
ever, the logarithm in this diagram comes from the region 
where the internal momentum carried by the transverse 
boson is sufficiently large such that one of the two vir- 
tual particle-hole pairs, either before or after the boson 
exchange, is far away from its mass-shell. This implies 
that the logarithm in the diagram (b) comes from two 
distant regions, and this exactly compensates the combi- 
natoric factor of two in the diagram (a). As a result, the 
logarithmic singularities from vertex and self-energy in- 
sertions to the polarization bubble cancel out, and there 
is no need to redo the calculations of logarithmical self- 
energy and vertex corrections. 



D. Renormalization group equations 

So far. wc have found singular logarithmic renormaliza- 
tions of the electron residue Z and the vertex due to the 
exchange of transverse bosons in one- and two-loop or- 
ders. The two-loop renormalization is exactly 1/2 of the 
square of the one- loop renormalization. One can attempt 
to go beyond perturbation theory and sum up the leading 
logarithmic corrections up to infinite order. The inter- 
play between one-loop and two-loop renormalizations is 
a good indication that the system is renormalizable such 
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that one can apply a renormalization group (RG) treat- 
ment. For this purpose, we introduce a running residue 
Z for the electrons and a running vertex F, defined at a 
certain energy scale D. Applying a standard procedure 
of integrating out high-energy degrees of freedom, one 
obtains RG flow equations for Z and F: 



g(l/^) _ 1 
91ogD 2kpa 
dV 1 



-ZV^ 



2t^3 



91ogD 2k pa 



Z^V 



(4.16) 

(4.17) 



The flow equation (|4.16p follows from the one-loop result 
for the self-energy, Eq. (|4.4p . after taking into account 
that a rcnormalizcd electron propagator and the rcnor- 
malized two vertices in Fig. |6l^a) yield an additional fac- 
tor ZT^ . Similarly, the RG equation (|4.17p follows from 
Eq. (j4.7p upon evaluating the diagram IHl^b) with renor- 
malized propagators and vertices. The flow starts at an 
energy scale Dq ^ Ep/{kpa) and is cut at either the dis- 
tance to the fermionic mass-shell or at the scale wpL of 
Eq. (|4.5p . The initial conditions for the RG equations 
are Z{Do) = 1 and F(Do) = 1- 

One can easily verify that the combination TZ is an 
invariant of the RG flow. 



TZ ^ 1. 



(4.18) 



This invariance is the manifestation of the cancellations 
between vertex and self-energy corrections which we ob- 
served in the explicit computations above. In partic- 
ular, the condition TZ = 1 ensures that the polariza- 
tions Hq, remain unchanged by the RG flow. Indeed, in 
the calculations with logarithmical accuracy, one can se- 
lect a cross-section with the smallest energies and add 
self-energy rcnormalizations and independent renormal- 
izations of the vertices on both sides of the selected cross- 
section (this is how we get an overall factor of 2 for the 
diagram [TUTb) in the previous subsection). This implies 
that the dressed diagram for the polarization bubble dif- 
fers from the bare one by {TZ)"^ and is obviously not 
renormalized. 

We emphasize that the relation TZ = 1 is valid to 
logarithmic accuracy, i.e., only the logarithmically singu- 
lar parts of the self-energy and vertex correction cancel. 
Non-logarithmic self-energy and vertex renormalization 
terms do not cancel, but arc small in our theory l/{akp). 
Exact cancelations due to Ward identities only arise in 
the limit of small momenta q/td — ^ 0. We are consid- 
ering, however, the opposite limit of small frequencies, 
uj/q — >■ 0, when there is no exact cancelation between 
self-energy and vertex diagrams. 

Using TZ = 1 the flow equation for the residue Z can 
be simplified to 



dZ 



1 



-Z. 



(4.19) 



i91og_D 2k pa 

This equation is easily solved and we finally obtain for 



the electron Green function 

G(k,.„)= ^[^"'"+Y"f-^-] (4.20) 

where the energy-dependent Z-factor is given by 

f f^y ifk|>..pL 

Zi^) = > ' < -n. (4-21) 

with 7/^ = l/{2kpa) and, Dq - Ep/{kpa). At the QCP, 



G(k,a;„i) oc 



r)-n± 



{iuj„i + S||(w,„)-efc)^"''^ 



(4.22) 



The form of the electron Green's function close to a ne- 
matic QCP, given by Eqs. (g^Hl), (H^, and (g^H), with 
the longitudinal self-energy I]|| defined in Eq. (|3.3p . is the 
main result of this work. Separating real and imaginary 
parts of G, we obtain Eqs. (|l.ip and (|1.2p presented in 
the Introduction. 

We emphasize that the Green's function contains the 
two energy scales uj^-^ ^ Ep/{{kpa)^X^) ex S,~^ and 
wpL ^ Ep/((kpa)'^X) ex £,^^ The self-energy Ey evolves 

at \oJm\ ^ WpL from a FL form at \ujm\ < Wp^ to a non-FL 

form at \uim\ > '^fl- ^^ addition, the Green's function 
develops an anomalous dimension r]± if the distance to 
the mass-shell exceeds the energy wpL. For kpa » 1 
which we assumed to hold, this anomalous dimension ri± 
is small. 



V. SUMMARY AND DISCUSSION 

We analyzed the electron Green's function in an 
isotropic metal in two spatial dimensions, close to a ne- 
matic QCP at which the Fermi sphere spontaneously 
develops a quadrupolar moment. The hallmark of this 
quantum phase transition is the presence of two criti- 
cal bosonic modes representing two quadrupolar polar- 
izations. These two modes are characterized by different 
dynamics. Whereas the polarization longitudinal to the 
quadrupolar momentum tensor is damped by particle- 
hole pairs in the metal and has the dynamical exponent 
zh = 3, the transverse mode remains undamped with 
z± = 2^ 

The self-energy correction due to the longitudinal fluc- 
tuations, E||, has been investigated beforei^"— '^"i"'^'^'^"' '— 
At one-loop order, its singular part depends only on fre- 
quency and acquires a strong non-FL form, E||(wm) ^ 

Ldm , for frequencies larger than |w„j| > cjp^ ~ C^^", 
where S, is the correlation length of the transition. In 
the present work, we focused on the modifications of the 
electron Green function due to the exchange of trans- 
verse bosons. We performed a systematic perturbative 
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analysis within an effective Eliashberg-type theory which 
operates with electron propagators already dressed by 
the one-loop longitudinal self-energy I]||, see Eq. (|3.3p . 
This effective theory has to be treated with extra care 
as it possesses spurious divergencies J^ For example, the 
transverse polarization at the one-loop order does not re- 
produce the original z^ = 2 dynamics of the transverse 
fluctuations. However, these divergencies are compen- 
sated by terms in the perturbative expansion which are 
formally of higher order in the number of loops and the 
fully renormalized bosonic propagator preserves z = 2 
dynamics. This has been first established in a pertur- 
bation analysis in Ref. [l^- We extended that analysis 
to higher orders and also included into consideration the 
curvature of the FS. We found that spurious divergen- 
cies indeed cancel out at all orders, and the curvature 
does not affect the cancellation. We found similar can- 
cellation of transverse renormalizations for the fcrmionic 
self-energy at the mass-shell. 

Our key result is the discovery that the exchange of a 
transverse fluctuation leads to a singular logarithmic cor- 
rection to the residue Z of the electron Green function 
already at one-loop order. The logarithm is cut by ei- 
ther the distance to the mass-shell or by the energy scale 
■^FL ^ C~^! whichever is larger. The apperance of the 
energy scale inversely proportional to ^ is an unexpected 
result because transverse bosons have z_l = 2, and scal- 
ing arguments suggests that a characteristic energy scale 
set by the transverse mode is Wp^ ^ C~^- This scale has 
manifestations in thermodynamics,— but we found that 
it is unimportant for the renormalization of Z arising 
from transverse boson exchange (see, however, the re- 
mark below on the renormalization of .J). The reason is 
that the logarithmic enhancement of Z comes from the 
region of the phase space where the transverse bosonic 
propagator is near its pole, i.e., J7,„ ~ vpq^ (aq), and the 
upper limit of the 1/q behavior of the integrand for Z is 
q ^ 1/a in which case r2,„ ~ vpq. The lower limit is £,~^, 
and the interplay between the two limits yields the z = 1 
scaling of ojfl- 

We also found, at one-loop order, a singular logarith- 
mic correction to the vertex F. We extended the calcula- 
tions to two-loop order and found log^ terms, both for the 
self-energy and the vertex. We then applied RG strategy 
and obtained flow equations for the running Z and F. We 
found that the flow equation satisfies ZT = 1 and also 
verified this result in explicit calculations to two-loop or- 
der. The condition ZT = 1 implies, in particular, that 
bosonic polarizations remain uncffected by logarithmical 
singularities although individual self-energy and vertex 



corrections to the polarization bubble are logarithmically 
singular. 

The solution of the RG equations yields our main re- 
sult: an electron Green function at the nematic QCP de- 
velops the anomalous dimension r]± = l/{2kpa), where 
kp is the Fermi momentum and a is a quadrupolar scat- 
tering length. The fully renormalized Green's function is 
given by Eq. (|l.ip . Away from QCP the anomalous di- 
mension persists if the distance to the mass-shell exceeds 
the energy scale wpL ~ ^~^- The anomalous dimension 
is detectable in, e.g., MDC ARPES measurements, see 

Eq. ini). 

There are additional logarithmically singular correc- 
tions in the theory that we neglected in our treatment 
for simplicity. First, it was recently argued by Metlit- 
ski and Sachdcv^ and Mross et ali^ that the longitudi- 
nal fluctuations also contribute a logarithmically singular 
correction to the fermionic propagator but only in three- 
loop order. It is an interesting open question as to how 
these additional logarithms exactly affect the renormal- 
ization group flow of the full theory and, more impor- 
tantly, whether they modify the zj^ = 2 dynamics of the 
transverse fluctuations or not. 

Second, in the mean-field theory the correlation length 
is given by .J ~ l/^yi -f gc,2, where gc,2 is the n ~ 2 
charge Landau parameter which approaches gc.2 ^^ ^1 
at the QCP. As shown in Ref. [ij, S, by itself acquires log- 
arithmic corrections when the self-interaction of nematic 
fluctuations is taken into account. In our present analy- 
sis, we also neglected these latter logarithmic corrections. 
We expect that the additional RG flow for the correla- 
tion length is decoupled from the flow of the fermionic 
residue, Z, so that only the expressions for the crossover 
energies wpL and Wp^ are affected. Nevertheless, due to 
this flow of ^ the frequency dependence of the residue of 
the Green's function should become sensitive to the so 
far elusive energy scale WpL ~ £,~^^ ■ 
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